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A (2 + l)-DIMENSIONAL GROWTH PROCESS WITH EXPLICIT 

STATIONARY MEASURES 

FABIO LUCIO TONINELLI 


Abstract. We introduce a class of (2 + l)-dimensional stochastic growth processes, 
that can be seen as irreversible random dynamics of discrete interfaces. “Irreversible” 
means that the interface has an average non-zero drift. Interface configurations cor¬ 
respond to height functions of dimer coverings of the infinite hexagonal or square 
lattice. The model can also be viewed as an interacting driven particle system and 
in the totally asymmetric case the dynamics corresponds to an infinite collection of 
mutually interacting Hammersley processes. 

When the dynamical asymmetry parameter {p — q) equals zero, the infinite-volume 
Gibbs measures (with given slope p) are stationary and reversible. When p ^ q, 
TTp are not reversible any more but, remarkably, they are still stationary. In such 
stationary states, we find that the average height function at any given point x grows 
linearly with time t with a non-zero speed: 'EQx{t) := ¥,{hx{t) — hx(0)) = V(p)t while 
the typical fluctuations of Qx{t) are smaller than any power of t as t —>■ oo. 

In the totally asymmetric case oip = Q,q = 1 and on the hexagonal lattice, the dy¬ 
namics coincides with the “anisotropic KPZ growth model” introduced by A. Borodin 
and P. L. Ferrari in HE]. For a suitably chosen, “integrable”, initial condition (that 
is very far from the stationary state), they were able to determine the hydrodynamic 
limit and a CLT for interface fluctuations on scale exploiting the fact that 

in that case certain space-time height correlations can be computed exactly. In the 
same setting they proved that, asymptotically for t —> oo, the local statistics of height 
fluctuations tends to that of a Gibbs state (which led to the prediction that Gibbs 
states should be stationary). 
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1. Introduction 

To motivate the object of our study, let us start with a well-known (l-l-l)-dimensional 
growth process. At all times t, the configuration is an integer-valued height function 
X G Z I— >■ hx{t) G Z with space increments hx — hx-i = ±1, see Fig. Local minima 
turn to local maxima with rate p (this corresponds to deposition of elementary squares) 
and local maxima to local minima with rate q (evaporation of elementary squares). If 
positive interface gradients are identified with “particles” and negative gradients with 
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Figure 1. The ASEP: squares are deposed (i.e. minima turn to max¬ 
ima) with rate p and evaporate (i.e. maxima turn to minima) with rate 
q. In the particle interpretation, particles jump to the right/left with 
rate q/p and cannot occupy the same site. 


“holes”, this process is equivalent to the one-dimensional Asymmetric Simple Exclusion 
process (ASEP). 

The study of this and similar stochastic growth processes in dimension (1 -|- 1) wit¬ 
nessed a spectacular progress recently, especially in relation with the so-called KPZ 
equation, cf. e.g. imiiiiiini for recent reviews. Some of the basic questions that were 
solved for certain models include the identification of the translation-invariant station¬ 
ary states (for ASEP, these are simply the combinations of Bernoulli measures for any 
intensity p G [0,1]), the determination of the dynamic scaling exponents characteris¬ 
ing the space-time correlation structure of height fluctuations, the study of the limit 
rescaled fluctuation process and its dependence on the type of initial condition. The 
same KPZ scaling relations appear also in the context of (1 -|- l)-dimensional directed 
polymers in random environment, last passage percolation and random matrix theory, 
just to mention a few instances [271 HI do]. 

On the other hand, for [d -|- l)-dimensional stochastic growth models, d ^ 2, the 
situation is much more rudimentary and mathematical results (see notably [Mill!) are 
rare. In this work we introduce a (2 -|- l)-dimensional stochastic growth process, for 
which we study the stationary measures and the corresponding large-time behavior of 
height fluctuations. The two-dimensional interfaces entering the definition of our pro¬ 
cess are discrete (i.e. heights are integer-valued) and are given by the height function 
associated to dimer coverings (perfect matchings) of either the infinite hexagonal or in¬ 
finite square lattice m- Height functions corresponding to dimer coverings of bipartite 
planar graphs, or to the associated tilings of the plane, are classical examples of discrete 
two-dimensional interfaces. Eor instance, dimer coverings of the hexagonal lattice (i.e. 
tilings of the plane by lozenges of three different orientations) correspond to discrete 
monotone surfaces obtained by stacking unit cubes, see Eigure “Monotone” means 
that if we let hx^y denote the height w.r.t. the horizontal plane of the vertical column 
of cubes with horizontal coordinates (x, y), then hx,y ^ max{hx+i^y, hx^y+i). In a sense, 
discrete monotone height functions are the most natural (2 -|- l)-dimensional analogue 
of the (1 -|- l)-dimensional height functions appearing in the one-dimensional ASEP. 


This work was partially supported by the Marie Curie lEF Action “DMCP- Dimers, Markov chains 
and Critical Phenomena”, grant agreement n. 621894. 
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Given a density vector p = {pi, P 2 , Ps) £ K+ with pi + P 2 + Ps = ^, there exists [22] 
a unique infinite-volume translation-invariant ergodic Gibbs measure vTp such that 

• the three types of lozenges have densities pi,i = 1, 2, 3 and 

• conditioned on the tiling configuration outside a finite region A of the plane, iTp 
describes a uniformly random tiling of A. 

The measures vTp have an explicit determinantal structure that will play a role in this 
work and that is recalled in Section [2.21 



Figure 2. A lozenge tiling of a portion of the plane, or equivalently a 
monotone stacking of unit cubes 

To model a growth process, we want to introduce a Markov evolution which is asym¬ 
metric or irreversible, in the sense that the interface has a net drift, proportional to an 
asymmetry parameter p — q. Moreover, as discussed in Section [l . 1 1 below, in order that 
its fluctuations can be at least heuristically described by a (2-|-l)-dimensional KPZ-type 
equation, the average interface speed should be a non-linear function of the interface 
slope. The most natural (2-|-l)-dimensional generalization of the ASEP described above 
(but which is not the one we will study here) would be the following. Let 

^'x,y • ~ hx^y-l) — hx^y ^0 ( 11 ) 

^x,y • — ^x,y Riax(/l3;_|_i y, hx^y-\-l) ^ 0, 

and observe that A+y (resp. A“y) is the maximal number of cubes we can add to (resp. 
remove from) column (x, y) while respecting the condition hx\y' ^ max{hx'+i^y', hx'^y'+i) 
for every (x',y'). For every column {x,y), we add a single cube with rate p if A+y > 0 
and remove a single cube with rate q if A~y > 0. In words, single elementary cubes 
are deposed (Fig. |^top) with rate p and removed (Fig. [^bottom) with rate q (com¬ 
pare with Fig. [^. We refer to this as the “single-flip dynamics”. If p = q there is no 
drift and the infinite-volume Gibbs measures vrp |22j are stationary and reversible. If 
instead p ^ q, the stationary states are not known, but they appear to be definitely 
very different from the equilibrium Gibbs measures [laiMiia]. This process has been 
studied numerically and one finds that typical interface fluctuations grow with time 
like t^, with (3 ~ 0.24... [T5l f^ . This is in sharp contrast with the ASEP, where the 
Bernoulli measures are stationary, irrespective of p being equal or different from q. In 
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the language of Section [Ml the two-dimensional single-flip growth process is believed 
to belong to the so-called isotropic (2 -|- l)-dimensional KPZ class when q ^ p- Unfor¬ 
tunately, the single-flip process is very hard to analyze mathematically and very little 
is known rigorously. 

In this work we study, instead of the single-flip dynamics, a different (2-|-l)-dimensional 
irreversible growth process, that we call “bead dynamics” for reasons that will be clear 
later (in the hexagonal lattice case, “beads” or “particles” correspond to horizontal 
lozenges as in Fig. j^. As discussed in Section iMl the bead dynamics belongs (in con¬ 
trast with the single-flip dynamics) to the so-called anisotropic (2-|-1)-dimensional KPZ 
class when q ^ p. Updates of the dynamics consist in adding or removing a random 


number ^ 1 of cubes at some column {x,y), in the following way (see Section 2.3 for a 
precise definition and Section 3.1 for the analogous construction on the square lattice). 
For every column {x,y), we assign 


hx^y -|- i for every i = 1,..., A+y (deposition of i 
hx,y — i for every i = 1,, A~y (removal of i cubes 


• rate p to the update hx^y 
cubes to column (x,y)); 

• rate q to the update hx,y - 
from column {x,y)). 

If p = q again there is no drift and the measures [22] are stationary and reversible. 
Somewhat surprisingly, vTp turns out to be stationary (but not reversible!) for any 
density vector p and for any value of p — q. This is the content of our hrst result. 
Theorem 2.4 The same then clearly holds also if we add to the generator of the bead 


dynamics the generator of another process w.r.t. which vTp is reversible. The measures 
TTp and their convex combinations are the only stationary measures that can be obtained 
as L —)• oo limits of stationary measures for the bead dynamics periodized on the torus 
of side L. In principle our result does not exclude the existence of other stationary 
measures that cannot be obtained this way; there might exist for instance analogs of 
the so-called “blocking measures” of one-dimensional asymmetric exclusion processes 

naEi. 

We emphasize that it is a non-trivial fact that equilibrium Gibbs measures should 
remain stationary in presence of dynamical irreversibility. As we mentioned above, this 
is false for instance for the single-flip dynamics. Typically, one expects that a Gibbs 
measure of a reversible dynamics remains stationary after introduction of a drift only 
when the reversible dynamics satisfies a so-called “gradient condition” |30ll2ni|2|. As we 
discuss in Section |4.1.1 , for the symmetric dynamics with p = q one can indeed identify 
a certain “gradient condition” that might help explain why Theorem 2.4 holds. 


It is important to emphasize that stationarity of the Gibbs measures means that, if 
the process is started from the distribution TTp, the law of interface gradients is time- 
invariant. However, overall the height function has a time-dependent random shift 
hxoii) — ^a:o(0) where, say, xq is the origin of the plane. On average hxo{t) — M(0) 
grows like {p — q)tV for some non-zero and slope-dependent V but the amplitude of its 
fluctuations cannot be deduced immediately from the stationary gradient measure iTp. 


Our second result, Theorem 3.1 says that the typical fluctuations of hx{t) — hx{0) grow 


slower than any power of t. Under a certain (technical) restriction on the interface slope, 
we can actually prove that fluctuations are at most of order -y/log t, which we believe 
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to be the optimal order of magnitude. Recall that, in sharp contrast, for the single- 
flip dynamics fluctuations were observed numerically P3ESI to grow like a non-trivial 
power of t. 

(stationarity of vTp). Checking stationarity is easy for the 


A word about Theorem 


2.4 


process periodized on the torus of size L, see Section]^ The extension to the infinite 
lattice is, however, non-trivial. One may expect that, when L is large, on local scales 
and for finite times the system does not feel the periodic boundary conditions and there¬ 
fore locally the dynamics on the torus and on the infinite lattice could be coupled with 
high probability. The situation is however more subtle: while on the torus the process 
is always well-defined, in the infinite systems one can easily construct initial configu¬ 
rations such that, for instance, beads (horizontal lozenges) escape instantaneously to 
inhnity. This is due to the fact that we allow for an unbounded amount of cubes to 
be deposed/removed at a time, since is not bounded. In order for the coupling to 
work, one needs to prove that for typical initial conditions and with high probability, 
the random variables A^,^ remain sufficiently tight in time during the out-of-equilibrium 


evolution. An important ingredient in overcoming these difficulties is the work [28] by 
Seppalainen on the one-dimensional Hammersley process dlESlIll. In fact, viewing 
beads as particles, the bead dynamics can be seen as a two-dimensional generalization of 
the Hammersley process, or more precisely an infinite collection of interacting Hammer¬ 
sley processes, see Fig. [^(a different two-dimensional generalization of the Hammersley 
process was introduced by Seppalainen in |29) : in that case a full hydrodynamic limit 
was obtained, but the stationary measures and the size of height fluctuations remain 
unknown). As a side remark, the single-flip dynamics can be instead visualized in a 
natural way as an inhnite collection of mutually interacting one-dimensional ASEPs, 
see caption of Fig. 

As we explain in some more detail in Section in the totally asymmetric case p = 
0, g = 1 and on the hexagonal lattice, the bead dynamics is the same as the interacting 
driven particle system introduced by A. Borodin and P. L. Ferrari in [H]^. In |3|, for 
a specific, deterministic initial condition, the hydrodynamic limit and the convergence 
of height fluctuations on scale \/logt to a Gaussian field were obtained. For such 
initial condition, the above-mentioned problem of proving that the dynamics is well- 
posed does not arise, simply because each bead has a deterministic, time-independent 
maximal position it can possibly reach, and therefore cannot escape to infinity. As we 
mention in Section on the basis of jH Prop. 3.2] it was natural to conjecture our 
Theorem 12.41 


I.l. Isotropic and anisotropic KPZ classes. In order to predict whether the fluc¬ 
tuations of a (2 -|- l)-dimensional growth process should be described by a KPZ-type 
equation, one should look at the Hessian of V, the average interface velocity consid¬ 
ered as a function of the interface slope. Indeed, the evolution of the fluctuations h 
in the stationary state of slope V4> should be governed on large space-time scales by a 
stochastic PDE of the type 

dth = i/A/i -|- Q{dxh, dyh) -|- white noise, (1-2) 

with V a diffusion coefficient and Q{-,) a quadratic form whose corresponding symmetric 
2x2 matrix is proportional to the Hessian of V at V0. (At present, it is not known 
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how to regularize such equation in order to make it mathematically well-defined, as was 
done recently for its one-dimensional analog Hi). 

The growth model is said to belong to the “anisotropic (2 -|- l)-dimensional KPZ 
class” when the two eigenvalues of the quadratic form Q have opposite sign, and to the 
“isotropic (2-|-l)-dimensional KPZ class” when they have the same sign. As discussed in 
[5] , the bead dynamics belongs to the anisotropic class (the eigenvalues can be computed 
explicitly from formula (3.5) below for V). 


Models in the anisotropic class are in a sense easier than those in the isotropic class. 
Indeed, in the former case it was predicted by Wolf |33j that the non-linearity Q is 
irrelevant as far as the large-time behavior of the interface roughness is concerned, i.e. 
the fluctuations of hx{t) — hx{0) should be of the same order i/logt as for the linear 
Edwards-Wilkinson equation m, where Q is set to zero. Theorem |3.1| and Eq. (|3.8[) 


confirm this prediction, for the bead model. Apart from the bead dynamics we study 
here, there are a few other (2 -|- l)-dimensional stochastic growth model models known 
to be in the anisotropic KPZ class, and all of them are exactly solvable in some sense. 
In this respect, let us mention the model introduced by Prahofer and Spohn in [26], for 
which height fluctuations are also known to grow like \J\ogt. See also |3l Sec. 3.3] for 
growth models in the same universality class; it would be interesting to see whether our 
result extend to these processes. 

The situation is very different for models in the isotropic KPZ class. In this case there 
are, to our knowledge, no exactly solvable models and only numerical simulations are 
available (see m for an overview). The non-linearity Q is expected to be relevant and 
to produce a non-trivial dynamical height fluctuation exponent. In particular, while 
neither the interface velocity V nor the stationary states of the (2 -|- l)-dimensional 
single-flip dynamics can be computed explicitly, the model is widely believed to belong 
to the isotropic KPZ class and, as mentioned above, the dynamical fluctuation exponent 
is numerically estimated to /3 ~ 0.24.. mm- 


2. Irreversible lozenge dynamics and stationarity of Gibbs states 


2.1. Configuration space. The Markov process we are interested in lives on the 
set of dimer coverings (perfect matchings) of the hexagonal lattice or equivalently 
the set of lozenge tilings of the whole plane. See Eigure[^ The “elementary moves” of 
the dynamics consist in rotating by an angle tt/3 three dimers around a hexagonal face, 
see Eigure [^ In this move, a horizontal dimer moves up or down a distance 1. The 
generic move of the dynamics (defined precisely in Section 2.3), that was described in 
the introduction as the deposition/removal of k cubes, can be seen as a concatenation 
of a random number A; ^ 1 of elementary moves in k adjacent hexagons in the same 
vertical column. We can therefore see each “horizontal dimer/lozenge” (we call them 
“beads” hereafteiQ) as attached to a “column” (an infinite vertical stack of hexagons): 
the bead can move up and down along the column but not change column. The set of 
possible bead positions can be identified with Z on, say, even columns and with Z-|- 1/2 


similar terminology was adopted in [B] for a model where bead positions take real values: such 
continuous model can be obtained from the dimer coverings of the hexagonal lattice in the limit where 
the density of horizontal dimers tends to zero, by suitably rescaling the lattice. 
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Figure 3. A portion of dimer covering of % and the corresponding 
lozenge picture. The two beads in column i are encircled and the vertices 
bo, Wo are marked. 



Figure 4. The two elementary moves. 


on odd columns. Note that beads of neighboring columns are interlaced: if on column 
t there are two beads at positions zi < Z 2 then necessarily in column — 1 there is a 
bead at a position zi < Z 3 < Z 2 and similarly for column i -\-l. 

Definition 2.1. For each dimer configuration a and bead b we let be the 

collection of available positions above it, i.e. positions that b can reach via a concate¬ 
nation of elementary moves that do not touch any other bead and do not violate the 
interlacing constraints. We define similarly Ifi = Ifi (a) as the collection of available 
positions below it. 


Remark 2.2. Given a finite or infinite subset A ofTi, we denote (t|a the dimer config¬ 
uration restricted to A and rj\A the configuration of beads restricted to A. If A = % we 
omit the index A. If (as will be the case in Theorem 2.4 below) every column contains 
at least one bead, a can be reconstructed by knowing just rj. In this case we will identify 
a dimer covering with a bead configuration. 
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2.2. Height function. On T-L, we take a coordinate frame where the axis ei forms a 
clockwise angle + 57 r /6 with the usual horizontal axis and the axis 62 an angle +' 7 r/ 6 , 
see Figure We also set 63 = —ei — 62 to be the vertical unit vector. 

Definition 2.3. Let LL* denote the dual graph of LL (it is a triangular lattice, whose 
vertices are vertices of lozenges). Vertices ofLL* are as usual identified with hexagonal 
faces ofLL. 

The height function h : %* ^ Tj is an integer-valued function, defined up to an 
arbitrary additive constant. When moving one step in the ei or 62 direction, the height 
increases by 1 when a dimer (or equivalently lozenge) is crossed and stays constant 
otherwise. 


Note that, with this convention, h corresponds to minus the height function with 
respect to the horizontal plane, and observe also that when moving one step in the 63 
direction, h decreases by 1 if no dimer is crossed and stays constant otherwise. 

Given p = {pi,p 2 ) £ with 0 < pi,p 2 < 1, and 0 < pi + /92 < 1 (we call p a non¬ 
extremal slope) there exists a unique translation-invariant ergodic Gibbs state TTp with 
slope p. This is a translation invariant probability law on the set of dimer coverings of 
LL, that satishes (cf. [211 Sec. 6 ]): 


• -Kp is ergodic with respect to translations by aei + he 2 , a, 6 G Z; 

• it satisfies the Dobrushin-Lanford-Ruelle equations: conditionally on the dimer 
configuration age outside a given finite subset Q C vTp is the uniform measure 
over all dimer coverings ag of Q compatible with age, i.e. such that age U ag is 
a dimer covering of LL] 

• it has slope p, i.e. 'Kpihx+si - K) = Pi, i = 1, 2. 

Note that pi is the density of south-east oriented lozenges, p 2 is the density of north-east 
lozenges and P 3 := 1 — pi — P 2 the density of horizontal lozenges. The non-extremality 
requirement on p means that all three types of lozenges have non-zero density. 

The measure Tip is in a sense completely known and has a determinantal representa¬ 
tion, that we recall here briefly (cf. in particular (2.3)), since it will be needed in the 
following. See [221 for further details. First of all, color sites of % white/black ac¬ 
cording to whether they are the left/right endpoint of a horizontal edge and let 'Hw,'hiB 
be the sub-lattice of white/black vertices. We denote wo,bo the black/white vertices 
indicated in Figure]^ and we let ^x,^x, with x = (xi,X 2 ) G T?, be the translation of 


Wo, bo by xiei -h ^ 262 . 

Take a triangle with angles 6 i = 7 rpj,i = 1,2,3 and let ki,i = 1,2,3 be the length 
of the side opposite to 0*. Define the Kasteleyn matrix K = {R"(b, w)}b 6 'Pg,we'Hvi/ 
as follows: If b,w are not nearest neighbors, then K{h,w) = 0. If they are nearest 
neighbors, then iF(b,w) = ki or k 2 or k^ according to whether the edge bw is oriented 
south-east, north-east or horizontal. 

Define also the matrix K~^ = {R"~^(w, b)}we'Hw,be'Ps 


1 1 f z dz dw 1 f z ^^'2 dz dw 

* (27ri)2 P{z,w) z w (27ri)2 /j, -\-k 2 W z w 

( 2 . 1 ) 


where the integral is taken over the two-dimensional unit torus T := {{z,w) G : 
\z\ = |r(;| = 1}. The long-distance behavior of K~^ is precisely known [22]: since the 
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polynomial P has two simple zeros on the torus, K ^ decays like the inverse of the 
distance so that in particular 


K Hwo,ba;)| ^ 


C{p) 

Xl\ + \X2\ + 1 


( 2 . 2 ) 


with C{p) <00 (this in general fails if p is extremal, e.g. if only one of the three dimer 
orientations has positive density). 

Given a set of (not necessarily horizontal) edges ei = (wi, bi),..., = (wfc,bfc) of 

Ti, the correlation function 7 rp{6e-^ .. (with 5^ the indicator function that there is a 
dimer at e) is given by 


7rp(<5ei ... (5eJ = ^fj-fir(bj,wi)^ det (K ^(wi,bj)) 


1 ^ ^ /c ' 


(2.3) 


Note, also in view of formula (2.1), that the r.h.s. of (2.3) is invariant if we multiply all 
ki by a common factor c, so that we may for instance fix the sum ki + k 2 + k^ to 1. 


2.3. Definition of the dynamics and stationarity of Gibbs states. The dynamics 
is informally defined as follows (cf. Fig. [^. To each column f and to each possible bead 
position z (horizontal edge oil-L) we associate two independent Poisson clocks of mean 
p G [0,1] and q G [0,1] respectively. We call them p-clocks and g-clocks, with obvious 
meaning. Clocks at different locations are independent. When a p-clock (resp. a o'¬ 
clock) at (f, z) rings, if z) is occupied by a bead we do nothing. Otherwise, we look 
at the highest (resp. lowest) bead (if any) on column i that is at position lower (resp. 
higher) than z: if it can be moved to z without violating the interlacing constraints 
then we do so, otherwise we do nothing. 

It is not obvious that the process is well defined on the infinite lattice. The danger 
is that beads could escape to -|-oo or to —00 in finite time (even in an arbitrarily small 
time). This may occur when spacings between beads in the initial configuration grow 
sufficiently fast at infinity. The problem is that the rate at which a bead moves, say, 
upward is p x |/^| and the average size of the jump is (|/^| + l)/2, and \I^\ is not 
bounded. 

Our first result (Theorem 2.4) is that the process is well defined for almost every 
initial condition sampled from TTp and that Tip is invariant. “Well-defined” means that 
the displacement of every bead with respect to its position at time zero is almost surely 
finite for every t ^ 0. In the symmetric case p = q, assuming that the process is well 
defined, invariance of the Gibbs measure is obvious because it is reversible. 


To precisely formulate the result, let us start by defining, given K = {Kp,Kq) G 
(M^)^, a cut-off process where p-clocks at distance more than Kp (resp. q'-clocks at 
distance more than Kg) from the origin of Ti are switched off. As long as Kp, Kg < 00 
there is no problem in defining the process on the whole Ti, since this is effectively 
a Markov jump process on a finite state space (once a particle is inside the ball of 
radius inax{Kq, Kp) it cannot leave it and therefore there is only a finite number of 
particles, determined by the initial condition, that can ever move). We call the 
configuration at time t, started from initial condition a. Given a column i, let Zt{i, n; K) 
be the position of its bead at time t, with zt{i, n; K) < zt{i, n -|- 1; K). The label n 
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Figure 5. A portion of the lattice with the allowed moves and the 
respective rates. Beads are drawn as black lozenges. When p = 1, g = 0 
or p = 0, g = 1 the process can be seen as an infinite set of discrete 
Hammersley processes, one per column, each interacting with the two 
neighboring processes. If instead we allowed particles to jump only by ±1 
with rates p/q, the process would be the single-flip dynamics described 
in the introduction and would correspond to an infinite collection of 
mutually interacting asymmetric simple exclusion processes (ASEP), one 
per column. 


is assigned in the initial condition and is attached to beads forever. For instance, one 
can assign the label (£, 0) to the lowest bead in ^ with non-negative vertical coordinate 
(in the initial condition). We assume hereafter that in each column there is a doubly 
infinite set of beads, i.e. the index n runs over all of Z. 

Two processes with different cut-offs K and K' can be coupled in the obvious way: 
their p-clocks (resp. q-clocks) are the same in the ball of radius va.m.{Kp,Kp) (resp. 
m.ui{Kq,K'q)). It is then easy to check that zt{£,n;K) is increasing w.r.t. Kp and 
decreasing w.r.t. Kg. We will then define 

zt{£,n) = lim lim zt{£,n;K) (2.4) 

Kq^OO Kp^OO 

to be the position of the {£, n)-th bead at time t for the process without cut-off. 

Assuming that zt{i,n) is finite for every call Xt the corresponding bead con¬ 

figuration and let be the law of the process {Xt)t ^ o started with initial distribution 
u (if u is concentrated at some a, then we write just P^). 


Theorem 2.4. For almost every initial condition sampled from Hp, with p 


p, ujooiv p a non¬ 
extremal slope, the limit (2.4) is almost surely finite for all {i,n) and t 0. More¬ 
over, TTp is invariant. More precisely, if f is a local bounded function of the dimer 
configuration one has for every t 0 




fifiXt)) = 17rp{da)E^f{Xt)=7rp{f). 


( 2 . 5 ) 
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Here, a function / is said to be local if it depends only on ua for some finite A. It is 
also possible to see (cf. Remark 7.7) that the limit \2A\ does not depend on the order 
how one takes the limits limA" ->.oo and limA" 


Theorem 2.4 
Section 


8.1 


IS 


J. q T ^ ±\. jy T ^ 

proven partly in Section ^ (existence of the dynamics) and partly in 


(invariance of vrp) 


Remark 2.5. With this result in hand, it is clear that one can construct many other 
driven processes that leave Tip invariant, simply adding to the generator of the bead 
dynamics another generator C with respect to which Tip is reversible (for instance, £ 
could be the generator of the single-flip dynamics with symmetric rates). 


It is a relatively standard fact to deduce from Theorem 2.4 that, if we start from vTp 
conditioned to have a bead say at the origin, then the law of the dimer configuration 
re-centered at the time-evolving position of this marked bead {tagged particle) is time- 
independent, see Section [8^ More precisely, fix a horizontal edge cq of 77. Given an 


initial condition a such that there is a bead at eg, call 4>t its vertical coordinate at time 
t (the horizontal coordinate does not change). Let also Xt := with Tx the 

vertical translation by a: G Z, be the dimer configuration viewed from the tagged bead 
and call Pj/ the law of the process {Xt)t ^ o started from some initial distribution v. 
Finally, let TTp be the Gibbs measure vTp conditioned on the event that there is a bead 
at eg. 


Proposition 2.6. The measure Up is invariant for the dynamics of the dimer configu¬ 
ration viewed from the tagged bead: for every bounded local function f and f ^ 0, 


frp{f{Xt))=np{f). 


( 2 . 6 ) 


3. Interface speed and fluctuations 

The stationary states vTp are characterized by an upward or downward flux of beads, 
according to whether p > q or p < q. The particle flux is directly related to the average 
height increase in the stationary state. While the height function was defined only up 
to an additive constant, one can define unambiguously the increase of the height at a 
face X from time 0 to t: Qx{t) ■= hx{t) — hx{0) equals the number of beads that cross 
the face x downward up to time t, minus the number of beads that cross it upward. 

For each horizontal bond e let b'^{e) (resp. b~{e)) be the lowest (resp. highest) bead 
in the column of e, at vertical position strictly higher (resp. strictly lower) than e. Also, 
call V{e,f) the collection of hexagons that b~{e) has to cross to reach position e and 
set V{e,'l) = 0 if this move is not possible (keeping the other beads where they are). 
Define V{e,\) similarly. 

The following result identifies the average height drift and shows that the fluctuations 
of Qx{t) in the stationary measure are smaller than any power of t: 

Theorem 3.1. For any face x, 

^np{Qx{t)) = t{q-p)J (3.1) 


with 


J = 7rp(|{e : x G I/(e,t)}|) > 0. 


(3.2) 
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For every 5 > 0, 

lim {Qx{t))\ ^ t^) = 0. 

i _irv-i ^ ^ 


(3.3) 


Note that only edges e in the same column as x and above it can contribute to J. 
The value of J is independent of x by translation invariance of Hp. The r.h.s. of (3.1) 
is linear in t because of stationarity of Tip and linear in (q — p) because the stationary 
state TTp does not depend on p, q. 

Theorem |3.1| is proven in Section 


It is not obvious to compute J explicitly in terms of the slope p, starting directly 
from the determinantal representation of the Gibbs states. In A. Borodin and P. 

L. Ferrari considered the dynamics for p = 0,q = 1 for a special, “integrable”, initial 
condition to, whose height function {ho{x))x£H* is deterministic and has non-constant 
slope (see Fig. 1.2 of [Ij: lozenges with a dot correspond to our south-east oriented 
lozenges, white squares to our north-east lozenges, while dark lozenges correspond to our 
beads). Let us emphasize that with such initial condition, each bead has a deterministic 
lowest position it can possibly reach on its column (this is related to the fact that in 
[U Fig. 1.1, 1.2] there is no dotted lozenge with coordinate n < 0), so that the well- 
posedness of the process poses no problem in that case. One of the results of |5 is a 
hydrodynamic limit, that in our notations we can formulate as follows: for every r] 
and r > 0 one has 


and h satisfies 


1 




drh. = V(d^h, drjh) 


1 sin(7rcl^h) sin(7r5^h) 
TT sin(7r(9r,h -|- cl.ch)) 


(3.4) 


(3.5) 


(this corresponds to formulas (1.9)-(1.11) in [3j, after after a suitable change of coordi¬ 
nates due to the fact that in BU the height is not taken with respect to the horizontal 
plane and a different reference frame than our ei,e 2 frame is used). From this, one can 
naturally guess that J in ( |3.2[ ) should be given by 

1 sin(7rpi)sin(7rp2) 

J =-(3.6) 

TT sm(7r(/9iP2)) 


Since pi,p 2 and pi + p 2 are in (0,1), the above expression is immediately seen to 
be positive. After a first version of this work was completed, Chhita and Ferrari [9| 
proved, through a smart combinatorial identity based on the determinantal structure 
of the Gibbs states, that indeed ( |3.6[ ) holds. 

By the way, Proposition 3.2 of [3] says that the law of local dimer observables around 
point (, [r]L\ ) at time tL tends as L —)■ oo to that of the same observables under 
the Gibbs state of slope p = ((9^h(^, r/, r), 9,,h(^, r?, r)). On the basis of this, it was 
natural to conjecture that our Theorem |2.4| holds. 


Referring to (3.3), we believe that the order of magnitude of the variance of Qx{t) 
is actually logt: this is indeed the result found by Borodin and Ferrari HE], in the 
particular case where p = 0, q = 1 and for the special initial condition ui mentioned 
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above. In this respect, our method allows indeed to refine estimate (3.3), under a (purely 
technical, we believe) condition on the slope p, to the following: 

Theorem 3.2. If the slope p satisfies 

\/%k^C{p) < 1 , 

with C{p) defined in ( A.25| ) and fci,/c 2 as in (2.1), we have for some c < oo 
limsupP^p(|( 53 ,(t) - ^ uy^logt) ^ 


t^OO 


(3.7) 


(3.8) 


For instance, if p = (1/3,1/3) (the three types of dimers have density 1/3, in which 


case ki,k 2 , k^ are all equal) one finds, evaluating numerically the integral in (A.25), that 


the l.h.s. of (3.7) is 0.896... < 1, so that (3.8) holds. By continuity, this remains true in 


a whole neighborhood of p = (1/3,1/3) while, again numerically, (3.7) does not seem 
to be satisfied in the whole set of non-extremal slopes p. 


Let us stress once more that we believe (3.8) to hold for every non-extremal p and 


to be of the optimal order w.r.t. t, while we do not attach any particular meaning to 


condition (3.7). 


Remark 3.3. It is possible to define alternatively the stationary drift as follows. Sample 
a from -kp and call as above fit the vertical coordinate of the tagged bead bo at time t. 
From Proposition fS.hl it is easy to deduce that the average of fit — fio is exactly linear 
in t, while from the definition of the process and the fact that \IfiJ has the same law 
as 17^1 (the Gibbs measures are invariant by reflection through the center of any given 
hexagonal face; this follows e.g. from uniqueness of-Kp given the slope) one sees 


1 


V := - I Ttp{da)E^(fit-fio) = (p-q)Ttp(\I()-fi{\I^J + l)/2). 


(3.9) 


^ I ru/ P\l bpl 

It is not hard to deduce from the stationarity of -Kp that 

IE,rp(hx(t) - ha;(0)) =-tpsu, (3.10) 

where we recall that p^ is the density of beads (the reason for the minus sign is that when 
a bead moves upward the height function decreases). Indeed, suppose for simplicity that 
p = l,q = Q. The l.h.s. of (3.10) equals minus the sum over the the edges e below x 


of the probability that there is a bead at e at time zero and that at time t is has moved 
at least n -|- 1 steps up, with n 0 the number of hexagonal faces between e and x. By 
translation invariance of-Kp, this equals 


-P3 


/ -kp{da)F„{fit - fio> n) = -ps / kp{da)E^{fit - fio) = -tpsv (3.11) 


n ^ 0 ‘ 


where we used positivity of fit — fio in the first equality and (3.9) in the second. 

3.1. Extension to dominos (perfect matchings of W?). Our result extends to per¬ 
fect matchings of T?, or equivalently domino tilings of the plane, cf. Fig. also in 
this case, one can define an asymmetric Markov dynamics (the height function has a 
non-zero drift) that leaves the Gibbs states invariant. We give only a sketchy descrip¬ 
tion of the generalization, omitting those details that are identical to the case of the 
honeycomb lattice. 
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ini 



Figure 6 . The correspondence between dimer covering (perfect match¬ 
ing) of and domino tiling of the plane. 


Since 7? is bipartite we can color its vertices black/white with the rule that each 
vertex has neighbors only of the opposite color. The height function h on the set of 
faces of 1 ? can be defined (modulo an arbitrary additive constant) as follows: for each 
X, y choose any nearest-neighbor path Cx^y from x to y and set 

hy-hx= ^ (Te((5e-l/4) (3.12) 


with the sum running over the edges crossed by the path, fje = ±1 according to whether 
Cx^y crosses e with the white vertex on the left/right and 5e the indicator function that 
there is a dimer at e. The definition is independent of the choice of path. 

The classification of translation-invariant ergodic Gibbs states is analogue to the 
honeycomb lattice case (actually the structure is the same for all planar, periodic, 
infinite bipartite graphs [ 22 ]): there exists an open polygon P C (for the lattice 

it is a square, while for P it is a triangle, as discussed in Section 2.2) such that for 
every p = (pi, P 2 ) £ P (non-extremal slope) there exists a unique translation-invariant 
ergodic Gibbs state vTp satisfying 


p{hx+ei — Pi, i — 1)2, (3.13) 

where the vectors e* are as in Figure The determinantal representation ( |2.3| ) still 
holds, with a different polynomial P{z,w) that however still has two simple zeros on 
the torus T. In order to define the irreversible dynamics that leaves the Gibbs states 


invariant, we have to find an analogue of the “columns” and “beads”. This is inspired 
by [211 [23| . The set of square faces of Z^ is sub-divided into infinite “columns” (indexed 
by £ G Z), i.e. diagonally oriented zig-zag paths, see Figure]^ Dimers that occupy an 
edge across a column are called “beads”. Each column is oriented along the positive ei 
direction, so it makes sense to say that a bead 61 in column i is above a bead 62 in the 
same column. 

Given columns £, f-l-1, call 1/ the set of vertices of 7? shared by the two columns and 
order the sites of 1/ according their ei coordinate. Then, a bead b on column i is said 
to be higher than a bead 6' on £ -|- 1 if the vertex of 6 on 1 / is higher than the vertex of 
b' on y^. With this definition, it is easy to see that beads satisfy the same interlacement 
property as on the honeycomb graph: given beads 61,62 on £, there exists 63 on ^ — 1 
and 64 on .^- 1 - 1 with 61 < 63 < 62 and 61 < 64 < 62. Also, like on the honeycomb lattice. 
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Figure 7. A portion of with vertices colored black/white. The 
colored graph is periodic, i.e. invariant w.r.t. translations by nici + 
71262,711,712 G N. The white and shaded zig-zag regions are the 
“columns”, indexed by i. Columns are oriented in the positive ei di¬ 
rection. The two thick bonds, transversal to columns, represent two 
beads. 


it is easy to see that if there is at least a bead in each column, then it is possible to 
reconstruct the whole dimer covering knowing only the bead positions. 

The dynamics is then defined as follows. Assign to any possible bead position, i.e. to 
each edge that is transversal to some column, two independent Poisson clocks of rates 
p and q, as before. All clocks are independent. When a p-clock (resp. g-clock) at edge 
e of column i rings, if there is a bead at e then do nothing. Otherwise, move the first 
bead below (resp. above) e in column £ to position e, provided this does not violate 
the interlacing constraints. Note that the dynamics is the same as on the honeycomb 
lattice, only the definition of “column” and “bead” being lattice-dependent. Observe 
also that each move can be seen as a concatenation of elementary moves on n adjacent 
faces along the same column, each elementary move consisting in the rotation by 7r/2 
of two dimers on the same face of Z^ (Fig. |^. In fact, the effect of an elementary move 


Figure 8. The two elementary moves of the domino dynamics 
is to shift a single bead one position up or down along its column. Note that, like in 
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the case of the hexagonal lattice, when a bead moves one step upward crossing a face 
/, the height function at / changes by — 1 . 

As in Section]^ given an edge e transversal to some column call b~{e) the highest 
bead in i, strictly lower than e and let l/(e,t) the collection of square faces of £ that 
b~{e) crosses when it is moved to e (with y(e,t) = 0 if the move is not allowed). Then: 

Theorem 3.4. The claim of Theorems \2.4\ and \3.1\ hold also for the bead dynamics on 
dimer coverings ofT?. 


With the exception of Section 4.2, in the rest of the work we will always consider the 
case of the hexagonal lattice. 


4. Dynamics on the torus 


4.1. Honeycomb lattice. We will let the torus denote the hexagonal graph TL, 
periodized (with period L) along directions 62,63 and we assume that L ^ 3. Note that 
now columns i along which beads move are L “circles” containing L hexagonal faces. 
We will say as before that a bead moves “upward” or “downward”, but what we mean 
is that it moves in the positive or negative 63 direction around the torus. 

Let Np be the set of configurations such that the height changes by [p 2 L\ (resp. 
[L/ 33 J — L) along any closed path winding once in the positive 62 (resp. 63 ) direction. 
On each column i there are [/O 3 LJ beads and bead positions on neighboring columns 
are again interlaced. We denote the uniform measure over Np . It is known that vr^ 
converges weakly to vTp, if the configuration space is equipped with the product topology 
|22] . Essentially, averages of bounded local functions converge. 


On Tl the process is defined similarly as in Section 2.3 for the infinite graph. For 
instance, when a p-clock at an edge e rings, one moves to e the hrst bead that is found 
when proceeding in the —63 direction from e along the same column, unless this move 
is forbidden by the interlacing constraint. The process is ergodic on Np , actually it is 
known that we can go from any configuration to any other by positive-rate elementary 
moves as in Fig. (see m Lemma 1] for details). 

Proposition 4.1. The measures TVp are stationary. 


It is actually easy to deduce, using ergodicity of the process in each of the sectors 
Np, that the only stationary measures are convex combinations of vr^ . 


Proof of Proposition f.l. Call the generator of the process. We want to check that 


■K^pC^ = 0 


(stationarity of tt^). One can decompose the generator as + C with 

involving only the up-jumps (related to the p-clocks) and C~’^ the down-jumps. It is 
sufficient to prove that TipC^’^ = 0, for the argument being the same. For every 
a G Np we have vr^ (cr) = \/\Np\. Given a G Njf let be the collection of a' G Njf 
that can be reached from cr by a single non-zero up-jump (not necessarily of length one) 
of a bead and let be the collection of a' G Np from which one can reach cj with a 

single non-zero up-jump of a bead. For every a' G we have , a) = p, while 
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= —p|r2o-| simply because the sum of row elements of the generator is zero. 
We see then 

[7r^£+’-^] (cj) = ^7r^(cj')£+’'^(CT',fT) = “ I^^D- (4-1) 

We want to see that ini- = |no-|. Note that |flo-| = Ylb while ini = Ylb \^b 

with the sum running over beads and being as in Definition 


2.1 


We will prove that 

:= |ni ^^1 — IDo-l is independent of a: as a consequence, it must be zero because the 


sum over cr of ( |4.1[ ) is zero. Assume that a' differs from cr only by a single elementary 
up-move of some bead b on some column i. Then, after the move the only beads b' 


are b itself and 

column £ + 1 that is “just above 6” (see Figure]^ and analogously for the others. It is 


that may have changed their values of are b itself and with the bead in 




bti 


♦ 


bi+i 


^rb7_ 




4- 


£-2 l-ll £+1 



a 


Figure 9. The configurations a and a' around b. In this example, once 
the bead b moves up by 1, IT decreases by 1, increases by 1 while 

0/) 1 I 1 


•^e+i 


I + and I^_ stay constant. 

°e+i °e-i 

clear that the contribution of b to Ho-/ — Ho- is +2: indeed, |I^| decreases by 1 and \I^\ 

increases by 1. Then look at column I — 1. One of the following two mutually exclusive 

cases occurs (Fig. 9): either |I^ | increases by 1 and |IT | stays constant or |I^ | 

|_f °e-i °i-i °i-i 

stays constant and j/T | decreases by 1. In both cases, the net variation of K^i — K„ 

from column £ — 1 is —1. The same holds for column I -|- 1 (since we are assuming 
L ^ 3, columns I ± 1 are distinct). Altogether, K^i — = 0. We have proved that 

Kfj is unchanged if we perform an elementary up-move. Given that the space state is 
connected, we proved that ATo- is constant (and therefore zero) on Np ■ □ 


^At the expense of being pedantic let us emphasize that, on the torus, the set of positions available 
“above” a bead means the set of positions reachable via moves in the -l-ea direction. 































18 


FABIO LUCIO TONINELLI 


The analog of Proposition 2.6 for the dynamics on the torus is the following. 


Proposition 4.2. Fix a horizontal edge eo on let a be a configuration such that 
there is a bead at eo and call fit be the vertical position of this bead at time t. Let vr^ be 
TTp conditioned to the event that there is a bead at cq. The law is stationary for the 
re-centered process 


Proof. The proof is very similar to that of Proposition 4.1 Call the part of the 


generator of the process involving only p-clocks. We have to show for every a 


7r^£+'^(a) = 0. (4.2) 

A symmetric argument then gives Trp£~’^{a) = 0. 

The measure is uniform among the \Np\ configurations with a bead at bq. We 
have equal —p times the number of configurations a' different from a that 

can be reached from a with a single move. The conhguration can change either because 
a bead different from bo (the bead that is at eo) moves, or because bo itself moves 
and then the dimer configuration has to be re-centered around the new tagged particle 
position. Note indeed that, when bo moves, necessarily the conhguration viewed from 
it changes, since the distance from the hrst bead above it decreases. The number of 
reachable conhgurations is then J2b^bo ~ '^b \^^\- Similarly, one sees that 


P 


7rg{a')C~^’^{a',a) = „ 



P 




Ei4- 


Then, the l.h.s. of (4.2) equals the r.h.s. of (4.1) (only with I/IA"^! replaced by l/|iV^|), 


that we know to be zero. 


□ 


4.1.1. A “gradient condition”. The bead dynamics on the torus has an trivial conserved 
quantity: the number of particles. There is however a less obvious one. For each of the 
L columns = 1,..., L dehne 




(4.3) 


with t he s um running over the beads of column 1. We have seen in the proof of Propo- 
that the “total charge” X = X^^'l is exactly zero. A simple computation 


sition 


4.1 


shows that, when p = q, the instantaneous drift of X^^'^ is 

hm ^E(AW(ai+ 5 ) - AW(a 4 )|a„s ^ t) = (Z^ - - zW)(ut)(4.4) 

5-s.O d 

with 


= -? E(I4,I(I4,I + 1)) + 5 E(l4,,l(l4fl + 1)). 


(4,5) 


This is a “gradient condition” |30j : the derivative of the charge at i is given by the 
divergence of a current, here Z^^'l — which is itself the gradient of a function Y 

of the configuration. 

As we mentioned in the introduction, conditions of this type are typically the key to 
guarantee that a reversible Gibbs measure remains invariant once an external driving 
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field that breaks reversibility is introduced, see e.g. [21 Sec. 2.5] and [20]. The unusual 
fact here (with respect to the more standard framework of e.g. the simple exclusion 
or zero range processes) is that the current associated to the local charge := 

~ \ does not seem to satisfy a gradient condition, while that of the non-local 

charge (integrated along the columns) does. Note that on the infinite lattice 
is not well-defined (it is just infinite). 


4.2. Square lattice. The hnite graph with periodic boundary conditions is defined 
like for the honeycomb lattice, except that the directions along which one periodizes 
are now ei,e 2 , see Fig. Note that each periodized column is a “circle” containing 2L 
square faces. The measure vr^ is defined as the uniform measure over dimer coverings 
of Ti such that the height changes by [TpiJ when winding once in the e* direction, and 
TTp (/) tends to vrp(/) as L —)• oo for every local observable / [ 22 ] . 

Like for the honeycomb lattice, one has 

Proposition 4.3. The measure is stationary. 

Proof. The only point where the proof differs w.r.t. the honeycomb lattice case is the 
way one shows that jno-j := ~ Ylb\^b\ 


as after (4.1). Recall that 


IS 


it is sufficient to show that, after any elementary move, the difference irio-l — jf^o- 
unchanged, whatever the initial configuration a is. 

When an elementary move is performed at a face / in column i, a bead b jumps from 
an edge e to P that has a common vertex with e. This common vertex belongs to either 
Yi or l£_i (recall that is the set of vertices common to columns P 1 ). Assume 
w.l.o.g. that the former is the case, as in Figure 10 and that P is higher than e in 


column i. After the move, |/()'| decreases by 1 and |/^ | increases by 1. On the other 
hand, it is clear that \I^\ is unchanged for beads b' on column ^ 1 , or on any other 

column except I — \. Therefore, we have to find a change -|-2 of jOo-j — jfla coming 
from column I — Call b~^, resp. b~, the first bead above (resp. below) b in column 
£—1, and call h' the bead “between” b'^ and b~ in column (. — 2. (The notion of ordering 
for beads in neighboring columns was introduced in Section [3.1[ ). Then, with reference 
to Fig. [T^ note that: 


• if b' is at or higher than edge 63 , then b'^ is at or higher than 65 and is the 

same, irrespectively of whether 6 is at e or e'. On the other hand, edges 67,68 
are accessible to b~ if b is at e' and are not if b~ is at 6 , so |/^| differs by 2 in 
the two cases. Altogether, when b is moved from e to e\ the contribution of h~ 
to the change of irio-l — is -|- 2 , as desired; 

• symmetrically, when b' is at or lower than 64 then b~ is at or lower than eg. 

When b is moved from e to e', does not change, while iv+l decreases by 2 , 

since 67 , 63 are not available positions any more. Again, we get a change -|-2 for 
irio-j — ^^ 1 , this time coming from 6 ^. 

• finally, suppose that b' is at ei or 62 . If 6 is at e then position 67 is available 
for 6 "*“ and 68 is not available for b~, while if b is at e' the opposite holds. As a 
consequence, both 6 + and b~ contribute -|-1 to the change of irio-j — jni- 


□ 




20 


FABIO LUCIO TONINELLI 



65 


e' . 







1 ^7 
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68 
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i- 

-2 

i- 

■ 1 



Figure 10 


Deducing stationarity of Hp on the infinite graph from stationarity of vr^ on the torus 
works exactly the same on % or Z^; for dehniteness, in Section]^ we will stick to the 
former case. 


5. The discrete Hammersley dynamics (DHD) 


On the way towards Theorem 2.4 let us switch for a moment to a one-dimensional 
interacting particle system known as Discrete Hammersley Dynamics (DHD) [H]. The 
configuration space of the DHD consists of particle configurations on Z (at most one 
particle per site). Each site of Z has an i.i.d. Poisson clock of rate 1. When a clock 
rings at a site x, if the site is occupied then nothing happens; otherwise, take the first 
particle to the right of x and move it to x. Note that each particle moves to the left 
with rate equal to the number n of empty sites before the next particle to the left, and 
the new position is uniform among the n sites. We call zt{n) the position of the 
particle (n G Z) at time t. Particles are labelled in the initial condition in such a way 
that ZQ{n) < zo(n + 1), with some arbitrary choice of whom to label 0 (for instance, it 
could be the first particle to the right of the origin). Labels do not change as particles 
move. 

The works [Il|28] consider instead the (continuous) Hammersley process [Tj, which is 
defined similarly as the DHD, except that particles live on M instead of Z: again, each 
particle moves to the left with rate equal to the available space before the next particle 
and the new position is chosen uniformly in the available interval. In |28j it is proven 
(among many other results): 

Theorem 5.1. //lim„_^_oo n~'^zo{n) = 0, then the dynamics is well defined at all times: 
the displacement of a particle with respect to the initial position is almost surely finite 
at all finite times. 


Theorem 5.1 extends immediately to the DHD |14) and is obtained with the help of 
a Harris-type graphical construction, that we recall here. To each site of Z associate 
an independent Poisson point process of density 1 on M'*': this is the set of times when 
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the clock at that site rings. Given a realization of all these i.i.d. Poisson processes and 
given 0 ^ s < t, —oo < a < b < oo, we can consider the set of all possible up-right paths 
in the rectangle (a, b] x (s, t], i.e. sequences (xi, ti),..., (x„, tn) of space-time points in 
the point process in the rectangle, with xi < • • • < and ti < ■ ■ ■ < tn- Note that 
inequalities are strict (for times this is not restrictive since with probability one there 
is at most one clock ringing at a given time). Let as in [U [28] L((a, s), (6, t)) be the 
maximal number of points of the Poisson processes on one such path. Let also 

r((a, s),t, k) = inf{/i ^ 0 : L((a, s), (a -|- h, t)) ^ k}. 

Then (this is given in [U [28| in the continuous Hammersley process where L and P are 
defined similarly, but the same holds true also for the DHD) for every t ^ 0 

zt{n) = ini {zo{j) +r{{zo{j),0),t,n - j)}. (5.1) 

j s£ n 

Note that the DHD has the following monotonicity property: 


Lemma 5.2 (Monotonicity for the DHD). If we take two initial conditions such that 
zoin) ^ ZQ(n) for every n and if we let them evolve using the same Poisson clocks, then 
the partial order is preserved at all later times. 


Proof. This is immediate from (5.1): if some zo{j) is changed to zo{j) — a, a G N, then 
r((zo(j), 0), t, n — j) increases at most by o. □ 


The representation (5.1) also allows to get an upper bound on the probability that 
the displacement of a particle is large. Indeed, if zt{n) — zo(n) ^ — k then there exists 
j ^ n such that 

mzo(j},0), (zoU) -k + (zo(n) - zo(j}),t}} ^ n-j. 

With a union bound, the probability (conditionally on the initial positions) that zt{n) — 
zoin) ^ — k is upper bounded by 

^ P(L{{0,0),{0,zo{n) - zo{j) - k,t)) ^ n-j), (5.2) 

j<n:zo(j)-zo{n) ^ -k 

where P denotes the expectation only with respect to the Poisson clocks. One ha^ 

P(L((0, 0), (h, t)) ^k)^ {th)’^/ik\f. (5.3) 

Indeed, there are h\/{k\{h—ky.) strictly increasing distinct sequences 0 < xi < ... Xk ^ h. 
Given one of these, the probability that there is an up-right path (xi,ti),..., {xk,tk) 
equals the probability that a Poisson random variable of average t equals at least k. On 
the other hand, if X is a Poisson variable of average t then for k ^ 1 

j-k 


p{x ^k)=Y, 


n ^ k 


n! 


E' 

m ^ 0 


-t 


^m+k 

{m -|- k)l ^ k \' 


(5.4) 


because (m + k)\ ^ klml. Then, (5.3) follows from 

h^ 


h\ 


1 




{h-k)\k\ k\ 


^see also | 28[ Lemma 4.1] that is given for the continuous Hammersley process 
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Let US call the law of the DHD started from an initial conhguration z = { 2 ^ 0 (^)Inez- 

t--j{zo{n) - zoU) - kr-i 


From (5.3) we have then 


^z{zt{n) - zo{n) ^ - k) ^ ^ 

j<n:zo{j)-zo{n) ^ -k 

This bound will be used in Section 0 


((n-i)!)2 


(5.5) 


6. The process started from tt^ is well-defined 

i.e. the bead displacement 


Here we prove “the first (and easier) half” of Theorem 2.4 


is finite for almost every initial condition sampled from -Kp. 

Proposition 6.1. Suppose that the initial configuration a is in the set 

Y = {a : for every i, lim n~‘^zo{£,n) = 0 = lim n~'^zo{i, —n)}. 

n^oo n^oo 

Then the process is well-defined at all times: for every {£,n), almost surely zt{i,n) — 
zo{£,n) as defined in (2.4) is finite for all t 0. 


( 6 . 1 ) 


Actually, when q = 0 (resp. when p = 0) the condition lim„^oo n~‘^zo{£, —n) = 0 
(resp. lim^^oo n~'^zo{£,n) = 0) is not necessary. Note also that i^piY) = 1 for any non¬ 
extremal slope. Indeed, ZQ{£,n) is just the sum of the first n inter-bead distances along 
column £. Since the measure vTp is ergodic for the action of Z^, n~^zo{£,n) converges 
TTp-almost surely to the finite limit l/ps. 


Proof of Proposition 6^, Fix some column £. We want to prove that, say, zt{£, n; K) — 
zo{£,n) is almost surely bounded away from minus infinity, uniformly in K and for 
every n. Take the DHD slowed down by a factor q (i.e. its clocks ring with rate q and 
not 1) with initial condition zo(n) = zo{£,n) for every n and couple the DHD and the 
bead dynamics by establishing that the g-clocks (within distance Kg from the origin) on 
column £ of the lozenge dynamics are the same as the corresponding clocks of the DHD 
(the DHD has no p-clock). Then, bead positions are dominated by those of the DHD, 
in the sense that zt{n) ^ zt{£, n; K) for all times and for all n. In fact, call Si,i 1 the 
ordered times when one of the finitely many clocks in column £ of the dynamics Xt-x 
ring. We have 2:o(‘) ^ zo{£,-;K) (actually with equality). At time sj" the inequality 
is still true, since the beads in £ have not moved while some DHD particles may have 
moved to the left. At time si, one of the following cases occurs: 

• a p-clock rings. Then, a bead might move upward and nothing happens for the 
DHD. We have in this case obviously 


^s+(-) ^ z+{£,■■, K). 


( 6 . 2 ) 


a g-clock rings at an edge e within distance Kg from the origin, but no bead can 


be moved to e without pushing other beads. Again (6.2) holds (for the DHD, a 
particle can move to the left). 

a g-clock rings at an edge e within distance Kg from the origin and the bead 
just above it, call {£,n) its label, can be moved to e. By assumption, for the 
DHD, the first particle at position greater or equal to e has index mfin. After 
the update, for the DHD one has particle m at e and for the lozenge process 
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one has bead {i,n) at e. All other particles/beads are unchanged. Clearly then 
(6.2) holds also in this case. 

The argument is then repeated inductively starting from time s^. 

Since, by Theorem 5.1, zt{n)—zo{n) > —oo almost surely, we conclude that zt{i, ri] K) — 
zo{i,n) is almost surely bounded away from —oo, uniformly in K. □ 

7. Large gaps and propagation of information 

Let Bji be the ball of radius R centered at the origin of R. 

Definition 7.1. Let A(R,t) be the largest integer n such that there exist horizontal 
edges ei G Br and e-i on the same column of ei, at distance n from it, such that at time 
t there is no bead between them. Also, let A{R, ^ t) = sup{A(i?, s), s ^ t}. 

We need a preliminary result, giving an upper bound on the probability of having a 
large gap among beads. We start from the case of the torus: 

Lemma 7.2. For G N there exists a constant C = Cq such that, for all T > 0, R ^ 1 
and L large enough. 


C, 


^‘f{da)¥^[A{R, ^ T) ^ C^logR] ^ T^. 


(7.1) 



7r^(A(i?,0) ^ C^logR) ^ C^R 


-C 

L 


(7.2) 


if C/^ is chosen sufficiently large. Using stationarity of tt^, this holds for every hxed 
t ^ 0. Then, 




^ Cc log R}d^ ^(T+1 )CqR 


(7.3) 


Let 


T = inf{t > 0 : A{R, t) > 2^^ log R} (7.4) 

and observe that, after time r, a clock has to ring in the ball logR before A(i?, t) 

becomes strictly smaller than log R (this is just a necessary condition; not every ring 
in Bji^C(;iogR decreases A{R,t)). Note that the realization of the Poisson clock rings 
at times t > r is independent of the process up to r (and of r itself). On the other 
hand, with probability u uniformly bounded away from zero, none of the 0{R?) clocks 
in logR rings in the time lag [r, r + 1/7?^]- In conclusion, 


El 

p 


cT+l 


l{A{R,t) ^ C^logR}dt 


T ^ T ] ujF? 


(7.5) 


Together with (7.3) we get that 

p,4r^r)^(r + i)Cci?2-Vn. 


(7.6) 
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We conclude by observing that {r ^ T} = {A(i?, ^ T) > 2C,^ logi?} and recalling that 
C can be chosen as large as wished. □ 

For the dynamics on T-L the same argument does not work since we do not know 
(yet) that vTp is stationary. A similar result however still holds, but the proof requires 
a comparison with the DHD we introduced above: 

Lemma 7.3. For any T < oo,( € N there exists a constant C = C{(,T) such that for 
R large 


17rp{da)F^[A{R, ^ T) ^ ClogR] if 


(7.7) 


A useful variant of Lemma 7.3 that we will use later (and whose proof follows almost 
exactly the same argument) is: 

Corollary 7.4. Fix a horizontal edge e and a time T > 0. For ^ 0 let 
the event that there exists a time s ^ T and horizontal edges on the same column as 
e, with e'^ at distance n"*" above e and e~ at distance n~ below it, such that at time s 
there is a bead at and no bead between them. There exists C = C{p,T) such that 

^.MT,n+,n-) ^ (7.8) 

For the proof of Lemma |7.3| we need the following preliminary result: 

Lemma 7.5. Recall that Tip is the Gibbs measure conditioned to have a bead at cq and 
that (ft — (fo is the displacement of the tagged bead at time t. Then, for every T > 0 
there exists a positive constant a = a{p,T) such that for every D ^ 0 

Fttp{3t ^T : \(j)t- (fol > F>) ^ aexp(-D/a). (7.9) 

Proof of Lemma |7.5[ To fix ideas let us prove that 

P^-p(3t ^ T : (ft - (fo ^ - D) ^ aexp(-L>/a). (7-10) 

We have seen in the proof of Proposition |6.1| that the downward displacement of a bead 
b is at all times stochastically smaller than the leftward displacement of a DHD particle 
(for the DHD with clocks of rate q) up to the same time, started from a configuration 
where the particles are at the same position as the beads in the column corresponding 
to b. Since the DHD particles move only to the left, the event ^ T :(ft — (fo^ — D} 
means that the DHD particle corresponding to b has moved more than D by the non- 
random time T. 

Call n the label of the tagged bead b in its column, initially at position zo{n) := eo, 
and go back to (5.5). Observe that if zo{n) — zo{n — r) ^ u then there are at most r 
beads in a set of u adjacent horizontal edges below zo{n). Using Lemma A.l| 
that, except with probability exponentially small in D, one has 

zo{n) — zo{n — r) < max [r, e{p)D] for every r ^ 1 


we see 


e(p) 


(7.11) 

for some positive £ depending only on the slope p. Then, from (5.5), on the event (7.11) 
FJMn) - Mn) < - D) < r(r/e(p)r 

r>e{p)D 

that decays super-exponentially in D. 


(r!) 




^ c(p,r)e-"'(^’^)^'°s^ (7.12) 


□ 
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Proof of Lemma 7.3. On the event A(R, ^ T) ^ Clogi? there exists a time s and 
a horizontal edge e G such that at time s there is no bead in the C log R horizontal 
edges immediately above or immediately below e. Assume w.l.o.g. that the former is 
the case. Let e+ (resp. e~) be the lowest horizontal edge above e (resp. the highest 
edge below e) where there is some bead b'^ (resp. b~) at time s. Call N ^ Clogi? the 
distance between e"*“ and e~. There are two possible cases: 

(i) at time zero bead b~^ is within distance A^/10 from e"*" and similarly b~ is within 
distance N/W from e~. This implies that at time zero there is no bead in a 
vertical interval of length N/2, centered on the face at distance N/2 above e. 
Since in the stationary measure the distance between neighboring beads has 
exponential tails (Lemma 


A.l) and N ^ Clogi?, this event has probability 


O (i?^ exp {—a{p) C log i?)) 


for some positive a depending only on the slope p, where the factor Rf comes 
from a union bound over all possible positions of e. Choosing C = Cq sufficiently 
large, we get a 0(i?“^) bound. 

(ii) At time zero, either b'^ is at distance n ^ iV/10 from e+, or b~ is at distance 
m ^ N/10 from e~. Say, to fix ideas, that the former is the case. This implies 
that at the (random) time s ^ T the bead 6+ has moved, say downward, a 
distance n ^ A^/10 with respect to the initial position. Thanks to Proposition 

the possible values of iV ^ Clogi? and over the the 0(i?^) possible positions of 
e gives a bound 0(i?“^) if C is chosen large enough. 

□ 


7.5, this has probability exponentially small in n. Summing over n ^ N/W, over 


As an application, we show that information does not propagate instantaneously 
through the system: if two initial conditions sampled from equilibrium differ only outside 
a ball of radius R, it is very unlikely that in a short time the discrepancy propagates to 
reach the center of the ball. It is useful to give a proof of this fact, since an extremely 
similar argument will provide the proof of Theorem |2.4[ For usual short-range systems 
one has a ballistic propagation bound: information does not travel more than a distance 
const X t in a time interval t (cf. for instance [251 Sec. 3.3]). The situation is more 
intricate here due to the presence of a-priori unbounded gaps among beads. 

Let the pair {a, a') G be distributed according to some law u such that a ^ -Kp, 
a' ~ TTp and a, a' coincide in Br. Couple the two processes by using the same Poisson 
clocks for both and call Fi, the law of the joint process (crt,cj(). Let 6e{t) G {0,1} 
(resp. (^e(t) G (0,1}) be the bead occupation variable at time t at a hxed horizontal 
edge e (say at the center of Br) for the process started from a (resp. a'). Let also 
Amax = max(A(i?, ^ T),A'{R, ^ T)), with A'{R, ^ T) referring to the process 
started from a'. 


Proposition 7.6. For every T < oo,C, G N there is a constant C such that 
lPi^(('5e(t))te[0,T] ^ (<^e(0)te[0,T]) ^ 


(7.13) 





26 


FABIO LUCIO TONINELLI 


Proof of Proposition \7.(^ We have from Lemma |7.3| 

Py((5e(i))te[o,r] ^ (<^e(^))te[o,T]) (7-14) 

^ IPi/ [(<^e(i))te[0,T] ^ {^'e{'t))t&[0,Tf ^raax ^ C'c +C(^/Ey‘. (7.15) 

On the event (5e(0)te[o,r] ^ ('^e(^))te[o,r] call ti the first time when 6e(ti) / S'^iti)- 
There are two possible cases: 

(i) (5e(t)~) = S'^{tf) = 1 and say 6e{ti) = 0 7 ^ = 1. In this case at ti a clock 

rings in the column 0 (the one of e) at a horizontal edge xi within distance 

^max from e and in configuration a' (but not in a ,-) a bead in a neighboring 

b 1 

column is preventing the bead at e to move to xi. At time ti there is therefore a 
horizontal edge ei in column ± 1 , with distance within Amax + l from xi (the +1 
is because it is on the neighboring column), where the bead occupation variable 
is different. 

(ii) Se{tf) = = 0 and say 6e{ti) = 0 7 ^ ^(.(ti) = 1. This means that at ti 

the clock at e rings (in this case we set xi := e) and that in configuration 
(but not in a'_) a particle in one of the two neighboring columns is preventing 

a certain bead (below e if the clock is a p-clock and above the origin e if it is 
a q-clock) to reach e. In particular, as in case (i), at time ti there is an edge 
ei in column ±1 within distance A^ax + 1 from xi, where the bead occupation 
variable is different. 

Call t 2 < h the first time s at which (iei('S) 7 ^ <5ei(s). On the event Amax ^ C(^logR, 
we have that t 2 > 0 because ei is in the ball Bfi where initial conditions coincide. 

We iterate the argument (cf. Fig. and as before we deduce that at t 2 there is 
an edge X 2 in the column of ei, within distance Amax from it, where a clock rings and 
an edge 62 in a column neighboring the one of ei, and at distance within Amax + 1 
from X 2 , where the bead occupation variable is different. The iteration stops when 
On is outside the ball Bn of radius R. Note that Xi,Xi+i are within vertical distance 
2{Amax + 1) ^ 3Amax and horizontal distance I from each other. 

Altogether, if (5e(t))te[o,r] ^ (<^e(^))te[o,r] then either Amax ^ C^logi?, or there 
exists: 


• a chain of sites xi,...,x„,, with Xj,Xj+i on neighboring columns, xi on the 
column £ = 0 and within distance C^logi? from e (the center of the ball Bn), 
\xi — Xj+i| ^ 3(7^ log ii and \xn\ ^ ii/2; 

• a sequence of times 0 ^ < tn-i < .. .ti ^ T such that either the p-clock or 

the g-clock at Xj rings at time t*. 

We get with a union bound 

IP.((5e(t))ie[0,T] # im)no,T]) ^ CJR^ + Y. 

n ^ R/(QCq logij) 


where Pn is the probability that a Poisson variable of average T{p + q) is at least n, 
while Nn is the number of all possible distinct chains xi,..., x^ of n sites with the above 
specified properties. Of course Nn ^ (C^logii)" for some constant while (5.4) gives 


p < p-»T.log(n/T(p+g))+n 
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X3 


Figure 11. An example of the iteration. We have that \xi — 
ei-i\ ^ Amax and je, - Xi-i\ ^ (A^ax + !)• When Xj+i = e* it means 
that we are in case (ii) above. The arrow follows the chain xi,X 2 , ■ ■ ■ 


The sum in (7.16) is o[R ^). 


□ 


Remark 7.7. Take a sampled from iTp and let he the coupled processes with the 

same Poisson clocks and the same initial condition, except that at has cutoff parameter 
K = {Kp,Kg) and a't has a different cutoff' parameter K' = {Kp,Kg). With the same 
ideas as for Proposition \ 7. h| it is possible to prove that 


TTp{da)F^{{6eit))te[o,T] ^ (5e(i))te[o,r]) = £{K,K') 


(7.17) 


with e{K,K') —)■ 0 when m.\n{Kp,K'Kq,K') —)■ oo. From this one can deduce that the 


order how the cutoffs in (2.4) are removed is irrelevant. 


8. Stationarity of Gibbs measures in the infinite graph 


We will prove Theorems 2.4 and 2.6 only for the honeycomb lattice. As for square 


lattice, once the result is proven on the torus (cf. Section 4.2), the extension to the 
inhnite system works exactly the same. 


.1. Proof of Theorem 2.4, Let us first of all prove (2.5) in the case where / = 


ni=i ^ei-, where ei,..., are horizontal edges [k G N) and 5e is the indicator function 
that there is a dimer at e. Choose R large enough so that all e* are in the ball Bn and 
say close to its center. Call vr (resp. vjff} the marginal of vTp (resp. Vp) on Bn (or, 
to be pedantic, on 77 n Bn) and let usp be sampled from vn and a'^^ from From 
convergence of to vip as L —>■ oo, we can choose L sufficiently large and a coupling of 
(ur, such that = a'^^ except with probability sn that tends to zero as R ^ oo. 

Let cr G kin and a' G (with the set of dimer coverings of T/,) be sampled 
as follows. The restrictions (<t_b^,cj^^) to Bn are sampled from {vnW^- Given the 
realization of {(yBnW'Bp)-' configuration (cj-^^^^, outside Bn are sampled 
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independently: from vrp(-|(TB^) and from tt^ (-Icj^^). We have therefore 


that a 


~ TT 


PJ 


a' ~ TT^ 


and they coincide in Br, except with probability er. 


Now couple the processes {cFt)t ^ O; ^ o started from cr, a' by establishing that the 
Poisson clocks in Br are the same for the two, while those outside Br are independent. 


Proceeding exactly like in the proof of Proposition 7.6 and using both Lemma 7.2 and 


7.3 to estimate the probability that A{R, T) ^ log R in any of the two processes, one 
hnds that, except with probability Er + 7?“^ = e'^, the bead occupation variables at all 
edges e*, i ^ A; for the two processes coincide up to time T. Therefore, for every t ^T, 


+ e'r — TTp (/) + e'r — TTpif) + e'r + Er 


( 8 . 1 ) 


where we used Proposition 4.1 (stationarity on the torus) in the second equality. Arbi¬ 


trariness of T < oo and of R proves (2.5) in the particular case / = n:=i (the larger 
T is, the larger we have to choose R and therefore L). 

When / is any bounded local function depending only on the conhguration of the 
horizontal dimers, it is always possible to write / as a finite linear combination of 
functions of the form OiLi the claim of the theorem holds also in this case. 

Finally, it remains to consider the case where / is a local function depending also on 
the configuration of non-horizontal edges. This requires a slightly different argument. 

Let us start with a simple observation, see Fig. 12 Let a, b be two horizontal edges 


b 



Figure 12. Given the horizontal dimers at a, b, c, d, all the dimer occu¬ 
pation variables at edges of hexagons si,..., is determined. In fact, 
all north-east oriented edges between a and c and between d and b are 
occupied by dimers, and the same is for south-east edges between c and 
b or between a and d. 


in the same column i and let si,... ,Sk be the hexagons of I included between a, b. If 
we know that the only beads in si,... Sfc are at a, b and if we also know the location of 
the two beads, one in column 7-1-1 and one in column 7—1, whose vertical coordinate is 
between that of a and of b, then we can reconstruct unambiguously the dimer occupation 
variables of all edges (not just horizontal ones) of hexagons si,... ,Sk- 
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Call S a finite collection of hexagons such that the union of their edges contains the 
support of /. Let Ak be the collection of hexagons that are at graph distance (on Ti*) 
at most K from S {S itself is a subset of Ai^). Let Ek be the event that, for every 
s € S, there are two beads in Ax, one below s and one above it. From the discussion 
above we know that, on the event Ex, the dimer conhguration on all the hexagons in 
S is uniquely identihed by the bead conhguration in K. Let 

giv) ■= ff(^lAx) = (8-2) 

which depends only on Let also Exit) the event that Ex is realized at every 

s ^ L We have 


E^^(/(w)) = lE^(i)) + omuF^^iExm) 

= E^^igi7jt)) + Oi\\f\\o,F^^iExm). 


(8.3) 

(8.4) 


From Corollary 7.4 we deduce easily that ex,t ■= ^-KpiExit)^) tends to zero as iL —oo, 
for every hxed ^Therefore, 


E^^ifiXt)) = E^^igirjt)) + 0(||/||ooeA,t) 

= T^pig) + 0(ll/l|ooeA,t) = TTpif) + 0(||/||ooeA,i) 


(8.5) 

( 8 . 6 ) 


where we used invariance of the Gibbs measure for functions of the bead conhguration 
in the second equality and 

TTpig) = TTpif) + OiWfWooSXfi) 

in the last (note sxfi ^ ^K,t)- We conclude by letting K ^ oo. 


.2. Proof of Proposition 2.6 This is very similar to the proof of Theorem |2.4t so 


we will be very sketchy. Given any R > 0 and e > 0 one can choose L sufhciently large 
so that there is a probability law for the random couple (cr, a') G x such that 
TTp and a, a' coincide, except with probability e, in the ball Br. This 


(T '^pi ^ 

is done like at the beginning of the proof of Theorem |2.4[ in fact, the total variation 
distance between the marginals on B^ of vip, vr^ tends to zero as L ^ oo (this is because 
the statement is true for the measures vTp, not conditioned to have a bead at the edge 
eo, and the probability to have a bead at cq is uniformly bounded away from zero). As 
in Theorem 2.4, given any o > 0, the coupled bead processes (at, a).) that use the same 


clocks in Br coincide up to time T in the ball Ba, except with probability e, provided 
R is larger than some iio(a, T). On the other hand, by comparing the displacement of 
a bead with that of a DHD particle, we see that if a is sufficiently large (depending 
only on T) the “tagged bead” stays within distance a/2 from its initial position up to 
time T, except with probability s. In conclusion, the processes Xt,X[ re-centered at 
the position of the tagged bead of crt,cj( coincide (except with probability 2e) up to 
time T in a ball of radius a/2 centered at the origin. Toge ther with the fact that the 
re-centered process X[ has law vr^ at all times (Proposition |4.2[) this implies the claim. 
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9. Speed and fluctuations 


Here we prove Theorem 3.1 and 3.2 about average speed and fluctuations of the 
growth process. 

Let A be the ixi box in % defined as the collection of hexagons obtained by translat¬ 
ing a fixed hexagonal face x (say, the one at the origin of H) by aei -|- 6 e 2 , 0 ^ a, 6 ^ 
Let 

Qkit) = - hx(0)) = ^ Qx{t). (9.1) 

Remark that 

QA{t + 5) = QA{t) - Y, -r)nA\+Y # A(e, I) n A| + RA,t,5 (9.2) 


with yf^ jy'S' the indicator that the p/g-clock at e rings once in the time interval [t, t+5], 
while the “error term” RA,t,5 includes the contribution to the change of Qa from the 
events where there are n ^ 2 edges ei,... ,6^ where clocks ring in the time interval 
[t, t + 5] and where, for every i ^ n, either |H(ej, |) n A| / 0 or \V{ei, t) n A| / 0. 


Proof of (3.1). We want to see that 

IE7r„[QA(t)] = {q-p)tfJ 


(9.3) 


with J defined in (3.2). By linearity we can assume i = 1, i.e. A = {x}. 

To see that R{x},t,6 can be neglected for (5 —)• 0 let b^{t) be the lowest/highest bead 
above/below x in the same column and let I{t) be the collection of horizontal edges 
included between b~{t) and b~^{t). Let also I{t,6) = U^g[i^i^ 5 ]/(s). Then, observe that 
the only clock rings that can contribute ±1 to Qxit + (5) — Qxit) necessarily occur in 
Then, 

^ ^ 2} (9.4) 

where A^ is a Poisson variable of average 6{p + <?)|L(t, b)\. Note that the law of |I(t, (5)| 


for the stationary process of law is independent of t and that, from Corollary 7.4 


the random variable |/(0, 5)| has exponential tails. Therefore, IE 7 rp|.R{x},t,( 5 l = 0{6‘^) and 
we see that 


^^7Tp[Qxit)] = -pTTpilie : X G y(e,t)}|) + gvrp(|{e : x G P(e,4,)}|) 

= {q-p)'^pi\{e : X G P(e,t)}|) 


(9.5) 

(9.6) 


where we used stationarity of TTp and in the last step its invariance by reflections through 
any hexagon. □ 


Proof of (3.3) and (3.8). We compute the variance of QA{t). We have (recall 
where again we can see that RA,t,5 ~ for 6 small with the same argument as above), 
letting for lightness of notation (•) := 


= 2{QA{t)Ki{at)) + {K2{at)) 


( 9 . 7 ) 









A (2 + l)-DIMENSIONAL GROWTH PROCESS 


31 


where 


Kn{a) = (-1)> |y(e, t) n Ar + g |y(e, |) n Aj' 


(9.8) 


and the sums run over all horizontal edges of %. We have then, recalling also (9.3), 
{K 2 {at)) = T^p{K 2 ) and dt{QK{t)) = 7 rp(A'i), 

-■= JtdQdt) - (QAiWf) = 2{(0A(t) - ((3A(*)))(A',(<r,) -t,(K,))) (9,9) 


+»a(A' 2 ). 


One has (see Appendix [A]) that 


and, for every 5 > 0 , 


sup£ "'vrpdATnl) ^Ci = Ci{n) < oo 

l 


- T^p{Ki)f] ^ C2{5) < oo. 


(9.10) 

(9.11) 

(9.12) 


Remark 9.1. It is likely that the variance of Ki is actually of order without any 
spurious correction. Indeed it is proven in [7] that, if f is a local dimer function and 
fx is f translated by x £ 1?, then {1/t) X)|a;| ^ eifx — T^pifx)] satisfies a CLT with finite 
variance. The problem with Ki is that | y (e, t) H A|, | l/(e, J,) n A| are not local functions. 
While in principle they are “almost-local” (the probability that they involve more than 
n dimers decays at least exponentially in n, see Lemma A.l]) , even proving the weaker 
(9.12) requires some non-trivial work. 


We have from (9.9), from stationarity and from (9.11), (9.12) 


^A/2(t) ^ ^/M2{t)^7:p[{Ki{at) - 7rp(Ki))2] + 2fCi 

^ + 2fCi 

from which it is then immediate to deduce that 

M2{t = £)^ 


Now we are ready to prove (3.3). Let xq be a face in A. Write 

^ - {Qx,m > \Qk{£) - (Qa(^))| ^ + 0 ( 1 ) 

where we used (9.15) to neglect the event that \Qk{(.) — {Qk{£))\ ^ . On the other 

hand 


(9.13) 

(9.14) 

(9.15) 

(9.16) 

(9.17) 


Qa{£) ~ {Qa{£)) — ~Ai + A 2 + A 3 

- ^ ^ [^x(0) hx^ifl') p{hx ^a:o)] 


(9.18) 

(9.19) 


isA 


+ '^[hx{£) - hxoii) - T^p{hx - hxo)] +f[Qxo(.^) - {Qxoi^))]- 

xGA 


(9.20) 
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(9.21) 


We have (see again Appendix [A|) 

TTp [A\] = Oi^Hogl) 
so that, using stationarity and Tchebyshev, 

\Ail\A2\^e\ogl, (9.22) 

with probability 1 — o(l). Finally, we note that if event ( |9.22| ) holds and at the same 
time \Qa{^) — (QA(f’))| ^ , one cannot have \Qxo(^) — {Qxo{^))\ ^ Eq. (3.3) is 

then proven (just let £ := t). 

Remark 9.2. If for a given slope p the condition ( |3.7[ ) is satisfied, then one can 
prove (cf. Remark A.2 below) that 7rp[(Ari — 7rp(A'i))^] = 0{i‘^log£), to be compared 
with ( 9.12[ ). Going back to ( |9.13[ ) one sees that dOl ) is then improved to M 2 {t = 
i) ^ C^Fmogi. Repeating the argument that starts with (9.16[), one sees immediately 


that (3.3) is improved into (3.8). 


□ 


Appendix A. Some equilibrium estimates 

Here we give upper bounds on the probability that, at equilibrium, there is a large 
gap between two consecutive beads in the same column. We use this information to 
deduce several useful equilibrium estimates. 

Let Jr be a set of r adjacent horizontal edges in the same vertical column of R and 
Nr the number of beads in Jr- 


Lemma A.l. Let p be a non-extremal slope. For every A > 0 and u > 0 there exists 
C = C{X,u, p) < oo such that, for every r G N, 

TTpdNr — P3r\ ^ ur) ^ C'exp(—Anr). (A.l) 


Recall that TTp{Nr) = psr. 


Proof. It is known (cf. |21[ Sec. 6.3]) that W is distributed like the sum of r independent 
but not identically distributed Bernoulli random variables Bi,i ^ r of parameter qi 


satisfying X)* Qi = rps and Y)i “ Qi) ~ (1/^^) logr as r oo 

. One has then 


Trp{Nr - psr ^ ur) = P{'^{Bi - qi) ^ ur] 

i ^ r 


(A.2) 

^ exp(— Aur) + (1 “ 


(A.3) 


i 


Since for every A > 0 there exists Ci = C'i(A) such that exp(x) ^ l + x + Cix^ for every 
X G [—A, A], we get for r ^ ro(p) large 


TTp{Nr - Psr ^ ur) ^ e ]J[1 + C 2 {X)qi{l - qi)] 

i 

^ ^ g-Awr+C' 4 (A)logr ^ C{X, u)e~^'^^^‘^. 


(A.4) 

(A.5) 


The claim is immediately extended to r ^ foip), possibly changing C to a new constant 
C{X,u,p). With a similar argument one estimates 'Kp{Nr — psr ^ — ur). □ 
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Proof of (9.11). Just note that |lJ(e,t)nA| = k implies that there are k hexagons just 
below e with no beads (or k + d of them, if e is at distance d from A), an event that 
has probability exponentially small in k (or k + d) thanks to Lemma A.l The average 
of Kji is then immediately seen to be of order □ 


Proof of (9.21). It is well known 
like the logarithm of \x — y\. Then 
estimate. 


that the variance of hx — hy under vTp grows 
a Cauchy-Schwarz inequality implies the desired 

□ 


Proof of (9.12). By Jensen’s inequality and symmetry it suffices to show that the vari¬ 
ance of 


fA = J2\V{eA)nA\ 


is 0(£^+'^). Write 

/a = /f + /f - /f := E + E n A| - E l^(c t) \ A| (A. 6 ) 

esA e^A esA 

and, again by Jensen, it is enough to estimate the variance of each of the three terms. 
(2) (3) 

This is easy for ^ and . Indeed, for instance if e is outside A and at distance de 
from it, then |IJ(e,t)nA| ^ n implies that there is a sequence of at least n + d^ adjacent 
hexagons starting from e, where no bead is present. This has probability exponentially 
small in de + n. As a consequence, if ei, 62 0 A then 


7rp[|I/(ei,t) n A| X |U(e 2 ,t) C A|] ^ cexp(-c'((iei -k dea)) 


(A.7) 
( 2 ) 


from which a bound on the second moment (and therefore the variance) of f)^ 

(3) 

easily follows. A similar argument works for /)( . 

The case of is much more subtle. Observe (cf. Fig. 13) that having |I/(e,t)| = 


n > 0 is equivalent to the following: the horizontal edge that is at distance n below 


6 is occupied by a dimer, and so are the n edges e\, 


e' 

) 


and e\,... 


i.e. 


1 


T(e,t)|=n - ^e- 

1=1 


(A. 8 ) 


/ \ 

A/ _\4 


Figure 13. The event F(e,t) = n for n = 3. 
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We have 

/a ^ - ^p(/a ~ (^-9) 

n>0 

with fn = ZlesA l|E(e,t)|=n- We use then Jensen’s inequality, {YjiUaif ^ if 

ti ^ 0) Yli^i = !> to got 

^p[(/a ^ n\p[{fn - TTp{fn)f] (A.IO) 

n>0 

(we chose tn = l/(C'n^) with C = Yln>o'^~‘^)- ft remains to estimate the variance of 
fn. Letting V{e;n) = l\v{e,t)\=n - 7rp(|F(e,t)| = n), write 

T^piifn -'^pifn)f] ^ ^ | 7 rp[y(e;n)F(e';n)]|. (A.ll) 

e,e'eA 


Since the event |y(e,t)| = n implies that there are n — 1 adjacent hexagons without 
beads under e, we have from Lemma A.l, for any A > 0, J > 0, 

| 7 rp[y(e;n)J^(e';n)]| ^ C7(A, J, 


Together with (A.ll) this gives 

T^p[{fn-T^p{fn)f]^C'{\,5,p) ^ | 7 rp[l/(e; n)F(e'; n) 

e,e'eA 


|l-5g-An 


(A. 12 ) 


where the constant A will be chosen later. 

Recall from (A.8) that l|y(e,t)|=n is a product of dimer indicator functions on a certain 
set of p = 2re + 1 (not all horizontal) edges of %. Call ei = (bi, wi),..., Cp = (bp, Wp) 
such edges and let Cp+i,..., e 2 p be the analogous edges corresponding to l|y(g/(of 
course Cj+p is just e* translated by e' — e). Now we use formula (2.3); 

7 rp[l/(e;n)R(e';n)] = (A.13) 


where det means that, since the variables V{e;n) are centered, when we expand the 
determinant in permutations a of { 1 ,..., 2 p} we have to keep only the permutations 
such that in the product there are N ^ 2 “special” terms of the type iL“^(wj, bo-(j)) 
with i ^ p and a{i) > p or viceversa (note N is always even). Thanks to (2.2), each of 
the N special terms is of order l/|e — e'\ for |e — e'\ large. We will consider therefore 
only the contribution of permutations such that N = 2 (those with N > 2 will give 
a sub-dominant correction when the sum over e, e' G A is performed; we skip details). 
W.l.o.g. we assume that the special terms are iL“^(wjj, bo-(jj)) and A'“^(wj 2 , bo-(j 2 )) with 
*iw(* 2 ) ^ P,'i' 2 ,cr{ii) > p (one has afterwards to sum over the O(p^) = O(n^) possible 
choices of ii,i 2 , cr{ii),a{i 2 ))- 

The contribution to det(iL“^(wi, bj))i ^ ij ^ 2 p from such permutations is 


^il,i2,cr{il),cr{i2) det(A" (Wj, bj))|]^ ^ j j p j 

X det(iL (wp bj)){p+l ^ i,j ^ 2 p,i^i 2 ,j^o-(ii)} 


(A.14) 

(A.15) 
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with £ = ±1 a sign that will play no role later. We claim that there exists C{p) < oo 
such that 

det(iir“^(wi,bj))i6/j6j ^ C{pY (A.16) 

for any r ^ 1 and sets I, J of cardinality r. If this is the case, from (|A.13) we have 




'f\n 


^[fci k2C{pYY^n‘^ 


o /|2 


(A.17) 


\e — e'l^ \e — e\ 

(recall that comes from the summation over the possible values of ii,i 2 , ( 7 (^ 2 )). 
Plugging into (A. 12) and choosing A sufficiently large we get 

1 


vrp[(/n - Mfn)?] ^ C"(A, 5 )e-("/ 2 )A ^ 


'eA 


|e - e'|2(i-<5) 




Using this estimate in (A. 10) we finally get 




(A.18) 


(A.19) 


as desired. The contribution from permutations with N > 2 gives instead 0{Y) since 
|e — e^|“2 is replaced by |e — e^\~^ that is summable over e’ ^T-L. 

It remains to prove ( |A.16 ). This is based on Gram-Hadamard type bounds (cf. for 
instance m App. A4]): if fi,gi,i ^ m are vectors in a Hilbert space and || • || is the 
norm induced by the scalar product (•,•), then 

\det{{fi,gj))ij ^rn\ ^ Yl \\fj\\\\9j\\- (A.20) 

j sg m 

The second observation (this trick is often used in constructive Quantum Field Theory, 
see again m App. A4]) is that one can rewrite ( |2.1[ ) as 

1 


K Hwx,ba;') = 


|•2^z j'2TT p—id{x2—X2)+i4>i^i—^l) 

I de / d(p— 

0 Jo 


{27ry Jo Jo Pie,cP) 

^ ^ i-^xiBx') 

y£l? 


(A.21) 

(A.22) 


where P{0, (p) = P(e*®, 2 : is the complex conjugate of a complex number 2 ; and 


^x(y) = 


1 


f27r 


f27r p-i0{x2-y2)+i(t>ixi-yi) 


(27r)2 


/ d9 d4>- 
lo Jo 


\HoA)\ 


r2TT 


Bx'iy) = 


(■2-k p-i9{x2-y2)+i<p{x\-yi). 


(272)2 


/O |P(M)| 3 /^ 


(A.23) 


(A.24) 


Finally one applies (A.16) together with the observation that ||Aa;(-)||, ||Ha;(-)|| are upper 
bounded by a constant. Indeed, 


\AY-)f = \\BY-)f = 


1 


f27r 1*277 

de I #- 


1 


(2vr)2yo Jo '|P(0,0)| 


=: Cip) 


which is finite since P has only simple poles on the torus. 


(A.25) 

□ 
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Remark A.2. For a given choice of p (and therefore of ki,k 2 ,k 3 ) it may happen that 
C{p) in (A.25) satisfies kik 2 C{p)‘^ < 1. In this case, from the first inequality in (A.17) 


together with (A.11) and (A.10) we see that 


^p[ifA^-T^pifA^)f]=0[ Yj 




:'eA 


e — e' 


/|2 


= 0{f log!). 


(A.26) 
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